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ABSTRACT 

We have derived the Galactic bulge initial mass function of the SWEEPS field in the mass range 
0.15 < M/Mq < 1.0, using deep photometry collected with the Advanced Camera for Surveys on the 
Hubble Space Telescope. Observations at several epochs, spread over 9 years, allowed us to separate 
the disk and bulge stars down to very faint magnitudes, F814IU « 26 mag, with a proper-motion 
accuracy better than 0.5 mas/yr (20 km/s). This allowed us to determine the initial mass function of 
the pure bulge component uncontaminated by disk stars for this low-reddening field in the Sagittarius 
window. In deriving the mass function, we took into account the presence of unresolved binaries, 
errors in photometry, distance modulus and reddening, as well as the metallicity dispersion and the 
uncertainties caused by adopting different theoretical color-temperature relations. We found that the 
Galactic bulge initial mass function can be fitted with two power laws with a break at M ~ 0.56 M 0 , 
the slope being steeper (a = — 2.41±0.50) for the higher masses, and shallower (a = —1.25±0.20) 
for the lower masses. In the high-mass range, our derived mass function agrees well with the mass 
function derived for other regions of the bulge. In the low-mass range however, our mass function 
is slightly shallower, which suggests that separating the disk and bulge components is particularly 
important in the low-mass range. The slope of the bulge mass function is also similar to the slope of 
the mass function derived for the disk in the high-mass regime, but the bulge mass function is slightly 
steeper in the low-mass regime. We used our new mass function to derive stellar mass-to-light values 
for the Galactic bulge and we obtained 2.1 < M/Lp&uw < 2.4 and 3.1 < M/Lpeoew < 3.6 according 
to different assumptions on the slope of the IMF for masses larger than 1M 0 . 

Subject headings: Galaxy: bulge — stars: evolution 


1. INTRODUCTION 

The knowledge of the stellar initial mass function 
(IMF) is a fundamental piece of information in many 
research areas of astrophysics. From a theoretical point 
of view, providing tight constraints on the IMF proper¬ 
ties in different stellar environments - both in the field 
and in star clusters - is mandatory to develop a co mplete 
and reliable theory of star formation (iMcKee fe Ostriken 
120071 and references therein). At the same time, from a 
phenomenological point of view, the IMF is a fundamen¬ 
tal property of stellar populations, and hence a crucial 
input in any study of galaxy formation and evolution. 
For instance, it represents an important ingredient in 
t he computations of Population Synthesis models (see 
lYazdekis et al.l 120151 and references therein), and hence 
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it affects our capability to extract the properties of stellar 
populations such as their luminosity evolution over time, 
the mass-to-light ratio, the total star formation rate at 
low and high redshifts, and so on. Therefore, it appears 
evident that to improve our knowledge of the IMF, or at 
least to have stronger observational constraints on this 
crucial ingredient, is of pivotal importance in many as¬ 
trophysics research fields. 

It is particularly important to analyze the properties 
of the IMF in various stellar environments such as the 
disk and the bulge of spiral galaxies in order to verify 
whether the well-known differences (in age and chemi¬ 
cal composition) in the stellar populations hosted by the 
distinct galact ic components have an impact on the IMF 
(|Zoccali et al.ll2003il . An additional reason that makes 
the study of the IMF in the bulge of spiral galaxies 
and elliptical galaxies important is due to the possibility 
that these spheroids could potentially contain the major¬ 
ity of t he stella r mass of the universe (see, for instance 
iFukueita et all (119981) . 

As we discuss, the IMF for the Galactic bulge is un¬ 
likely to be very different from the present-day mass func¬ 
tion (PDMF) below the main-sequence turn-off (MSTO), 
since most of the star-formation in the Galactic bulg e 
happened within about 2 Gyr (IClarkson et al.l 120081) , 
with no evidence of star formation after that. So we 
will refer to the observed PDMF of the Galactic bulge 
as the IMF in the mass range below the MSTO, which 
occurs at «1.0 Mq for a stellar popul ation with sol ar 
meta llicity and an age of t = 11 Gyr (Calamid a et al.l 
12014 hereafter Paper I). 
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In spite of the huge improvements achieved in the 
observational facilities, there is not yet any chance to 
directly measure the IMF of spheroids outside of our 
Galaxy. As a consequence, the measurement of the IMF 
in the Galactic bulge is a fundamental benchmark (or ref¬ 
erence point) for any analysis devoted to in vest ig ate this 
prope rty in extra-galactic spheroids dCalchi Novati et all 

l 2008 lh 

Concerning the Galactic bulge, the two most recent de- 
terminatio ns of the IMF in our s pheroid have been per- 
for med bvlHoltzman et a.l~l (|1998l hereinafter H098) and 
bv iZoccali et al.1 (|200flL hereafter ZOOO), by taking ad¬ 
vantage of the exquisite observational capabilities of the 
Hubble Space Telescope ( HST ). In particular, the analy¬ 
sis performed by ZOOO pushed a step forward the knowl¬ 
edge of the bulge IMF thanks to the use of the Near- 
Infrared Camera and Multi-Object Spectrometer (NIC- 
MOS) available at that time: the derived mass function 
still represents the deepest measured to date and extends 
to ~ 0.15 Mq. They found a power-law slope for the 
IMF equal to a = —1.33 (when a Salpeter IMF would 
have a = —2.35, where dN/dM = C x M“.), with some 
hint for a possible change of the power slope - a ~ —2 
at ~ 0.5 Mq. ZOOO also found that the derived bulge 
IMF is steeper than that measured for the Galactic disk 
(|Reid &; Gizislll997t iGould et al.lll997f). In t his co ntext, 
it is also worth noting that iDutton et all (|2013f) have 
recently used strong lensing and gas kinematics to inves¬ 
tigate the existence of possible differences in the prop¬ 
erties of the IMF between the disk and the bulge in a 
sample of spiral galaxies within the Sloan WFC Edge-on 
Late-type Lens Survey dTreu et alJl20lU SWELLS). As 
a result they found a significant difference between the 
bulge IMF and that of the disk, the former being more 
consistent with a Salpeter IMF, and the latter being more 
consistent with a Chabrier-like IMF. 

On the basis of this evidence, it appears quite appropri¬ 
ate to analyze the properties of the Galactic bulge IMF 
in different fields of view and using more u pdated obser¬ 
vatio nal datasets. In a previous paper (1 Calamida et al.1 
12014 hereinafter Paper I), we have taken advantage of 
the availability of a huge photometric dataset for the 
low-reddening Sagittarius window in the Galactic bulge 
collected with the Advanced Camera for Survey (ACS) 
on board HST , to obtain the first unambiguous detection 
of the white dwarf cooling sequence of the Galactic bulge. 
In this manuscript, we use the same data to perform a 
thorough analysis of the IMF in this field of the bulge 
in order to provide additional, independent insights on 
the bulge IMF, thus supplementing the results of previ¬ 
ous analyses. In this investigation we explore a larger 
and denser field compared to what was previously ob¬ 
served by H098 in the Baade’s Window and by ZOOO 
in a field at (/ = 0.277 °,,b = —6.167°,). Most impor¬ 
tantly, for the first time we estimate the Galactic bulge 
IMF based on a clean sample of bulge stars thanks to 
the very accurate proper motions (down to FS1AW ss 
26 mag) that we were able to measure. Furthermore, we 
use a statistical approach to apply a correction for the 
presence of unresolved binaries. We note that the slope 
of the very low-mass range of the IMF is fundamental to 
estimate the mass budget of a stellar population, since a 
major fraction of the stellar mass is included in this range 
and low-mass stars have been hypothesized to contain 


a sig nificant fractio n of the total mass in the universe 
(|Fukuedta et al.lfl998fl . It is even more important in the 
case of the Galactic bulge since this component might 
include ps 20% (1.8 x 1Q 10 Mm) of the m ass of the Galaxy 
(jSofue et, alJl2009t iPortail et al.ll2015l ). 

The structure of the paper is as follows: in §2 we 
discuss the observations and data reduction in detail, 
while in §3 we describe how we selected a clean sample 
of bulge stars. In §4 we present the theoretical mass- 
luminosity relations we adopted to convert the luminos¬ 
ity functions in mass functions, while §5 deals with the 
different systematics that affect the estimate of the ini¬ 
tial mass function. In §6 we compare the derived IMF 
for the bulge with the disk mass function, and in §7 we 
derive a minimum value for the stellar mass-to-light ra¬ 
tio of the Galactic bulge. §8 deals with the gravitation 
microlensing events predicted by the bulge IMF derived 
in this work, and the conclusions are presented in §9. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the Sagittarius Window Eclipsing Ex¬ 
trasolar Planet Search (SWEEPS) field (l = 1.25°, b = 
—2?65) in the Galactic bulge in 2004 and again in 2011, 
2012 and 2013 with HST , using the Wide-Ficld Channel 
of ACS (proposals GO-9750, GO-12586, GO-13057, PI: 
Sahu). The SWEEPS field covers « 3.'3 x 313 in a region 
of relatively l ow extinct ion in the bulge ( E(B — V) < 
0.6 mag; lOosterhoff fc Ponsenlll968 1. The 2004 observa¬ 
tions were taken in the F606W (wide V) and F814W 
(wide J) fi lters over t he cou rse of one week (for more 
details see iSahu et all 1200611 . The new data were col¬ 
lected between October 2011 and October 2013, with 
a ~ 2-week cadence, for a total of 60 F606W- and 61 
F814IU-band images. The two datasets, the 2004 and 
the 2011-2012-2013 (hereafter 2011-13), were reduced 
separately by using a software that performs simultane¬ 
ous point-spread function (PSF) photometry on all the 
images. The choice to reduce the two datasets separately 
is due to the high relative proper motions of the disk and 
bulge stars in this field, caused by the Galactic rotation: 
the disk star relative proper motions (PMs) peak at pi ss 
4 mas/yr, with a dispersion of « 3 mas/yr, whereas the 
bulge motions are centered at m ~ 0 mas/yr, with a dis¬ 
persion of ~ 3 mas/yr (see Paper I). This means that 
a substantial fraction (~ 30%) of stars would move by 
more than half a pixel (25 mas) in 9 years. 

We calibrated the instrumental photometry to the 
Vegamag system by adopting the 2004 photometric zero- 
points, and we obtained a catalog of « 340,000 stars 
for the 2004 and for the 2011-2013 datasets. The left 
panel of Fig. 1 shows the F814W, (F606TV — -F814W) 
Color-Magnitude Diagrams (CMD) for all the observed 
MS stars in the 2011-2013 dataset. 

The right panel shows the sample completeness as a 
function of the F814IU magnitude. Details on how the 
completeness was derived are given in i f2.ll This figure 
shows that the completeness is ~ 50% at F814W ~ 25.5 
mag. The completeness of the E606VF magnitude is ~ 
50% at F606W ~ 28 mag. The 2004 dataset has a very 
similar completeness, reaching ~ 50% at F81AW ~ 25.3 
mag and F606W ~ 28 mag, respectively. 

In order to obtain a clean bulge MS sample to derive 
the IMF, we estimated the PMs of the stars in this field 
by combining the astrometry and the photometry of the 
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2004 and the 2011-2013 datasets. By comparing the po¬ 
sitions of stars in the two epochs we estimated PMs for 
« 200,000 stars down to F81AW « 25.5 mag. 

2.1. Artificial star tests 

To properly characterize the completeness of the mea¬ 
sured PMs, the photometric errors and the errors due 
to the reduction and selection techniques adopted, we 
performed several artificial star (AS) tests. We created a 
catalog of ~ 200,000 artificial MS stars, with magnitudes 
and colors estimated by adopting a ridge line following 
the MS. We then produced a second artificial star cata¬ 
log, by using the same input magnitudes and colors, but 
applying a PM to each star. We assumed the bulge PM 
distribution as measured in Paper I, with the distribu¬ 
tion centered at /q ~ 0 mas/yr, and a dispersion of ~ 3 
mas/yr. Artificial stars were added and recovered one by 
one on every image of the two datasets by using the same 
reduction procedures adopted earlier. In this way the 
level of crowding is not affected. In order to estimate the 
magnitude and color dispersion of the MS due to photo¬ 
metric errors and data reduction systematics, we selected 
recovered artificial stars with A Mag = (M agi — M ag 0 ) < 
0.5 mag, and d = ^{X a - A,;) 2 + (Y„ - l)) 2 < 0.75 
pixel, where the quantities with subscript i represent the 
input, and o represent the output, in both datasets, end¬ 
ing up with a sample of 146,225 stars. We applied this 
selection because stars which were not recovered in a cir¬ 
cle of radius 0.75 pixel can be safely considered not found. 
The left panel of Fig. 2 shows the selected recovered ar¬ 
tificial stars for the 2011-2013 dataset (red dots) plotted 
in the F8UW, (F606W - F814W) CMD; the observed 
stars are plotted as well as grey dots. The right panel 
shows the recovered color spread of the MS as a function 
of the F814W magnitude. The comparison of the artifi¬ 
cial (red dots) and observed (grey dots) CMDs indicates 
that we are not able to reproduce the entire color spread 
of the MS by assuming the presence of a single stellar 
population of solar metallicity and age t = 11 Gyr. A 
metallicity spread of more than 1 dex is present in the 
SWEEPS bulge field based on medium-resolution spec¬ 
tra of MS turn-off, sub-giant and red-giant branch stars 
collected with FLAMES at the Very Large Telescope (see 
Paper I for more details). The metallicty spread can fur¬ 
ther broaden the MS, and differential reddening, depth 
effects as well as binaries might also play a role. It is 
worth mentioning that most stars in the color and mag¬ 
nitude ranges 2.0< F606W — F814W < 2.5 and 18.5 
< F81AW <21.5 mag belong to the (closer) disk popu¬ 
lation. 

In order to estimate the completeness of the measured 
PMs, we matched the two recovered sets of artificial stars 
and compared the output with the input PMs in the 
direction of both X and Y axes as a function of the 
two magnitudes. Fig. 3 shows this comparison for the X 
(top panel) and the Y (bottom) axes versus the F814W 
magnitude. Only stars with A Mag < 0.5 mag and d < 
0.75 pixel are shown. This plot shows that the dispersion 
of the recovered PMs increases at fainter magnitudes as 
expected and the accuracy of the measured PMs is better 
than 0.1 mas/yr (« 4 km/s at the distance of the Galactic 
bulge) at magnitudes brighter than F81AW < 23. At 
F81AW ~ 25 mag where the completeness is > 50% for 


both datasets, the recovered PM scatter is ~ 0.25 mas/yr 
(« 10 km/s) within 3 er uncertainties. This precision will 
allow us to separate bulge from disk stars down to very 
faint magnitudes and to characterize the Galactic bulge 
mass function down to the very low-mass (VLM) range, 
M < 0.5 Mq. 

3. A CLEAN BULGE MAIN-SEQUENCE SAMPLE 

We adopted the measured PMs to select a sample of 
main-sequence (MS) stars devoid of disk-star contami¬ 
nation from the 2011-2013 dataset. PMs are projected 
along the Galactic coordinates and we considered stars 
with /i; < — 2masyr _1 to belong to the bulge, following 
the criteria adopted in Paper I. This selection allowed us 
to keep ss 30% of bulge members while the residual con¬ 
tamination of the sample by disk stars is < 1%. We ended 
up with a sample of 67,765 bulge MS stars. Note that the 
total contamination by disk stars in the S WEEPS fiel d is 
» 10% , as shown in the previous work of lClarkson et all 
(2008). Fig. 4 shows the F606W (left panel) and the 
F81AW (right) versus (.F814IV — F606W) CMDs of the 
selected bulge MS stars. The magnitude range covered 
by MS stars is different when observing with the -F606W 
or the -F814W filter, decreasing from ~ 8.5 to 7 magni¬ 
tudes. This happens because very low-mass MS stars are 
cooler and thus more luminous at longer wavelengths. 

The CMDs of Fig. 4 show that the color spread of the 
MS did not substantially decrease compared to the CMD 
of Fig. 1, confirming that most of the color dispersion is 
due to the spread in metallicity, the presence of some 
amount of differential reddening, depth effects and bina¬ 
ries. 

Fig. 5 shows the observed PM-cleaned bulge MS lumi¬ 
nosity function (dashed line) based on the F81AW mag¬ 
nitude for the stars plotted in the CMDs of Fig. 4. The 
completeness measured from the AS test is used to cor¬ 
rect the number of observed stars per magnitude bin and 
the corrected luminosity function is over-plotted in the 
same figure with a solid line. We applied the complete¬ 
ness correction by binning on the observed magnitudes; 
in this way we take into account the uncertainties due 
to the photometric errors moving the stars among the 
magnitude bins. The F814W-band corrected luminos¬ 
ity function of Fig. 5 extends from just below the bulge 
MSTO at F81AW ~ 19 mag down to F81AW ~ 26 mag, 
where the completeness level is ~ 30%. A similar lu¬ 
minosity function is obtained by adopting the -F606W 
magnitude. 

4. THE MASS-LUMINOSITY RELATION 

In Paper I we used the BaSTI 0 dPietrinferni et al.l 
12004 1200611 stellar-evolution database to fit isochrones 
to the bulge CMD. In order to extend the BaSTI 
isochroncfo to the range of very-low-mass stars 
(M < O.5M 0 ) we computed very-low-mass (VLM) stel¬ 
lar models for exactly the same chemical composition 
of the BaST I ones, by adopting the same physical in¬ 
puts used in 'Cassisi et al.1 (1200(1 hereafter CA00). We 

9 http://albione.oa-teramo.inaf.it/ 

10 In their standard format the minimum initial mass in the 
BaSTI isochrones is equal to O.oMq . The BaSTI isochrones ex¬ 
tended in the VLM star regime are available at the BaSTI URL 
repository. 
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Fig. 1.— Left: F814VF, (F606W — F814W) CMD of MS stars in the SWEEPS field based on the 2011-2013 dataset. Error bars are 
also labelled. Right: completeness of MS stars as a function of the FS14IF magnitude. The horizontal lines in both panels represent the 
F814W magnitude at which the completeness is 50% and the vertical line in the right panel represents the 50% completeness level. 
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Fig. 2.— Left: F814VF, (F606VF — F814W) CMD of recovered artificial stars for the 2011-2013 dataset (red points). Observed stars are 
marked with grey points. Stars are selected in magnitude, A Mag < 0.5 mag, and in position, Ad < 0.75 pixel. Right: (F606VF — F814VF) 
photometric scatter as a function of the F814VF-band magnitude as estimated from the artificial star test. 
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Fig. 3.— Comparison of Input and Output proper motions in 
the X (top panel) and Y (bottom) axes of stars recovered in the 
AS test as a function of the .F814VF magnitude. The 3cr limit is 
indicated by the overplotted red dots. Dashed and black dotted 
lines mark a dispersion of 0.1 and 0.5 mas/yr, respectively. 


note that, the accuracy and reliability of the BaSTI 
models and their extension to the VLM stellar regime 
have been extensively tested by comparing them with 
observed CMDs and mass-luminosity (M-L) datasets for 
both field and cluster stars. As a result, a very good level 
of agreement has b een obtained with the various observa- 


tional constraints (Cassisi et al. 200C 

: Bedin et al.. 20091: 

Cassisi ct al. 2009) Cassisi & Salaris 

20lit Cassisi et al. 

2014|). Since the VLM stellar models have been com- 


puted by using a different physical framework compared 
to the models of more massive stars in the BaSTI li¬ 


brary (see iPietrinferni et ah (20041 ) and CA00 for more 
details on this issue) one can expect that, in the stellar 
mass regime corresponding to the transition between the 
BaSTI and the VLM stellar models occurring at about 
~ 0.6 Mq, some small mismatch in surface luminosity 


and effective temperature at a given mass is possible. 
Since in retrieving the IMF one has to rely on the first 
derivative of the theoretical M-L relation, it is important 
to eliminate any such discontinuity in the M-L relation 
(jKroupa fe Toutl[l997D . To this aim, we devoted a huge 
effort - which included computing additional stellar mod¬ 
els using both the physical inputs adopted for the BaSTI 
library and that used by CA00 - in order to match the 
two model datasets at the stellar mass with (almost) the 
same luminosity and effective temperature. 

The evolutionary predictions were transformed from 
the theoretical to the observational plane by adopting 
the color-T e ff relations and bol ometric correction scale 
for the ACS filters provided bv lHauschildt et all (|19990 
for T e ff < 10,000 K, whil e at l arge r T P ff we adopted the 
relations published bv lBedin et al.1 ( 2005 1. 

Fig. 6 shows selected scaled-solar isochroneiO for 
the same age, t — 11 Gyr, and different metallicities, 
Z= 0.008, 0.0198, 0.03, plotted in the F814IV versus 
log{M/M & ) plane. We selected models with this age 
and abundances based on the fit of the bulge CMD per¬ 
formed in Paper I (see Fig. 2) and on the spectroscopic 
metallicity distribution for this field. In the same plot 
a solar metallicity isochrone but for an age of 8 Gyr is 
also shown (blue solid line). As expected, in the explored 
age and stellar mass range, the M-L is completely unaf¬ 
fected by an age change. In order to check the impact 
on the adopted M-L relation related to the use of a dif¬ 
ferent bolometric correction scale, we also plotted the 11 
Gyr, solar metallicity isochrone transferred in the obser¬ 
vational plane by using the standard Johnson bo lometric 
corrections provided bv IPietrinferni et all ( 2004 ) and the 
transformatio ns from the Jo hnson to the HST photomet¬ 
ric system bv lSirianni et all (120051 red solid). 

The five mass- luminosity relations all show a slight 
change of the slope around log(M/Mo) fs -0.3 (M ss 
0.5 Mq). This inflection point is due to the molecular 
Hydrogen recombination occurring at a mass equal to 
ss 0.5M 0 ; the formation of the H 2 molecule changes the 
value of the adiabatic gradient and, hence, the stellar 
structure thermal stratification (see Cassisi et al. 2000 
and references therein). 

Fig. 6 also shows the impact of using various metal¬ 
licities or ages for the selected M-L relation. As dis¬ 
cussed, for old ages, t > 8 Gyr, suitable for the Galactic 
bulge population under scrutiny, the exact value of the 
selected age is quite irrelevant. On the other hand, the 
change in the mass derived (at a fixed magnitude) using 
two different mass-luminosity relations corresponding to 
Z = 0.008 (which is the most metal-poor chemical com¬ 
position we selected) and Z = 0.03 (which is our most 
metal-rich composition) is only about a; 0.04-0.08 Mq 
in the high-mass range (M > 0.5 Mq), and « 0.02-0.04 
Mq in the lower-mass range. The spectroscopic metal- 

11 Our referee correctly pointed out that bulg e stars appear to 
be a-enhanced up to about solar metallicity lIZoccali et al.l 120081 ; 
IGonzalez et al.1120111 ; IJohnson et aljfeOlll . 12013) . However, we de¬ 
cided in present work to adopt scaled-solar models due to the lack 
of suitable alpha-enhanced VLM star sequences in a wide metal¬ 
licity range. This notwithstanding, we note that all the compar¬ 
isons performed in present paper are performed at constant global 
metallicity (and not at constant [Fe/H] and it is well known that 
a-enhanced stellar models are nicely mimicked by scaled-solar one 
with the same global metallicity (see, e.g. Pietrinferni et al. 2006 
and references therein). 
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licity distribution we derived for the SWEEPS field, as 
discussed in §3.1 and Paper I, spans a range of metal- 
licity from [M/H\ ~ -0.8 to ~0.6, i.e. more than 1 
dex. However, the distribution shows three peaks at 
[M/H] ~ -0.4, 0.0 and 0.3 and most of the stars, ~ 
85%, are included in the range -0.5 < [M/H] < 0.5. 
We can then safely assume the aforementioned metal- 
licity values, Z = 0.008andZ = 0.03, corresponding to 
the more metal-poor and the more metal-rich peaks of 
the distribution, to test the effect of metallicity on the 
mass estimate. However, we also tested the effect of fur¬ 
ther decreasing the metallicity of the adopted models, 
by using an isochrone for Z = 0.002 and the same age, 
t = 11 Gyr, to convert luminosities into masses. In this 
case, the mass estimate changes by ~ 17% in the entire 
mass range, when going from the more metal-rich model, 
Z = 0.03, to the more metal-poor, Z = 0.002. For a 
small fraction of stars in our field, less than ~ 10%, the 
mass estimate will have a ~ 5% larger uncertainty. 

The impact of using a different bolometric correction 
scale for transferring the models from the theoretical to 
the observational plane in the derived masses is smaller 
and of the order of ~ 0.005 Mq in the entire mass regime. 

5. THE GALACTIC BULGE INITIAL MASS 
FUNCTION 

The mass-luminosity relation we obtained for MS stars 
by using BaSTI isochrones is only the first step towards 
determining the IMF of the Galactic bulge. Uncertain¬ 
ties due to the assumed distance and reddening, presence 
of differential reddening, metallicity dispersion, depth ef¬ 
fects, and the presence of binaries need to be taken into 
account. 

Following iSahu et all (|2006T l and Paper I, we fitted the 
bulg e CMP using a distance modulus of /r 0 = 14.45 mag 
(ISahu et al.li200d ) and a mean reddening of E(B — V) = 
0.5 mag and a standard reddening law. Extinction coeffi- 
cients for the WFC fi lters are estimated by applying the 
iCardelli et all (|1989f ) reddening relations and by adopt¬ 
ing a standard reddening law, Ry = Ay/E(B—V) = 3.1, 
finding Afqoqw = 0.922 x Ay, Apguw — 0.55 x Ay, 
and E(F606W - F814W) = 1.14 x E(B - V). It is 
worth ment ioning that i f we us e the reddening value es¬ 
timated by iNataf et~ahl (120131) for the SWEEPS field, 
E(V — I) = 0.79, and their extinction curve, Ry = 
A V /E(B -V) = 2.5, we obtain E(B - V) = 0.47, in 
good agreement with the value we assumed. 

We used the T814W-band luminosity function to 
probe the bulge IMF since MS stars are brighter at redder 
colors and so the photometry in this filter is more com¬ 
plete and accurate than in the F606W filter for the same 
mass (see the CMDs in Fig. 4). We converted observed 
magnitudes into masses using the mass-luminosity rela¬ 
tion for solar metallicity, Z = 0.0198, and for an age of t 
= 11 Gyr, transformed by using the color-T e ff relations 
by lHauschildt et all (11999H . As we showed in the pre¬ 
vious section, age does not significantly affect the mass- 
luminosity relation for t > 8 Gyr, and observational evi¬ 
dence shows that most bulge star s in our field are older 
than 8 Gyr ([Clarkson et al.ll2008 . Paper I). 

In order to estimate the effect of dispersion in metal¬ 
licity, we computed the difference in the masses de¬ 
rived by adopting three different metallicities: solar 
(Z, = 0.00198), metal-rich (Z = 0.03), and metal-poor 


(Z = 0.008). For magnitudes in the range 18.0 < 
T814IU < 26, this changes the inferred masses by 0.02 
to 0.08 Mq, resulting in an uncertainty of ~ 8% in mass. 

We also varied the assumed distance modulus by 0.2 
mag, from 14.35 to 14.55 mag, corresponding to a depth 
of ~ 1 Kpc, and the extinction, E{B — V), from 0.45 
up to 0.55 mag. Both the distance and reddening un¬ 
certainty affects the derived masses by ~ 0.01 Mq over 
the entire mass range, i.e. 1-5%. Similarly, adopting dif¬ 
ferent color-Teff relations changes the derived massed by 
less than 2%. 

By summing in quadrature the uncertainties related to 
the parameters of the bulge CMD as fitted to our data, 
including metallicity, distance, reddening, and color-T e ff 
relations, we obtain a final systematic uncertainty on the 
mass estimate for each star. This uncertainty varies de¬ 
pending on the inferred mass and is carried through the 
remainder of the analysis. 

5.1. The effect of unresolved binaries 

Unresolved binaries, i.e. the expected presence of equal 
or lower-mass binary companions for many of the main 
sequence stars we observe, are likely to affect the inferred 
IMF of the bulge esp ecially at lower ma sses, such as 
M < 0 .5 Mq (see, e.g.. lKroupa et all (|1991l) and iKroupal 

<HHD). 

The availability of photometry in two different filters 
for a large fraction of our stars potentially allows us to 
correct for the effect of unresolved binaries, as binary 
systems will be somewhat redder than single stars of the 
same apparent brightness. However, the photometry, in 
particular at the faint end, is not sufficiently accurate for 
a direct identification of individual binary systems; our 
correction must therefore be probabilistic. 

Both the fraction of binaries and the distribution of 
mass ratios for the Galactic bulge are not well con¬ 
strained. However, in Paper I we showed that there 
is evidence for a substantial fraction of He-core white 
dwarfs in the bulge based on the color dispersion of the 
cooling sequence and the comparison between star counts 
and predicted evolutionary lifetimes. According to stan¬ 
dard stellar evolution models, He-core white dwarfs can 
only be produced in a Hubble time by stars experiencing 
extreme mass-loss events, such as in compact binaries. 
Indeed, in Paper I we reported our finding of two dwarf 
novae in outburst and five candidate cataclysmic vari¬ 
ables in the same field, both of which are characteristic 
of a population of binaries. Our evidence at the time sug¬ 
gested that the Galactic bulge has a fraction of binaries 
of larger than 30%. 

For the present analysis, we assume that the distribu¬ 
tion of mass ratios of binar y stars in th e bulge follows 
the distribution derived by IDuauennov fc MavoH (119911 
hereafter DM) for a sample of 164 F- and G-dwarf stars in 
the solar neighborhood. The distribution is a log-normal 
and it is described by the functional form: 


in the interval [0,1], where q = M 2 /M\, [i = 0.23 and a q 
= 0.42 and C = 10,900 for our sample of bulge stars. 

We also repeated the experiment by assuming a flat 
mass-ratio distribution similar to the distribution found 


m = Cel 
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Fig. 4.— Left: PM-cleaned bulge MS F606VF, F606W — F814VF CMD; note that 70% of the bulge stars were rejected because of the 
PM selection. Error bars are also labelled. Right: Same stars plotted on the F814W, F606W — F814PF CMD. 



F814W [MAG] log (M/M 0 ) 

Fig. 5.— .F814VF-band observed luminosity function for the Fig. 6.— Theoretical mass-luminosity relations for different 

Galactic PM-cleaned bulge MS stars (dashed line) and the lumi- metallicities, ages, and color-temperature relations, 
nosity function corrected for completeness (solid line). 


by iRaghavan et al.l (|2010t ) based on data for 454 F- to 
K-dwarf stars within 25 pc of the Sun. 
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Fig. 7. — IMF of the Galactic bulge. The two power laws that fit 
the IMF are over plotted, for a slope a = -2.41±0.50 (dotted line) 
and a = 1.25±0.19 (dashed), together with a log-normal function 
with M c = 0.25 and a = 0.50 (green solid). The Salpeter mass 
function (blue dashed-dotted line) and the Chabrier log-normal 
function (red dashed double dotted) are also shown. Error bars 
are displayed. 


In a simplified Bayesian approach, we use the fraction 
of binaries and the distribution of mass ratios from Equa¬ 
tion (1) as a prior for the presence and mass of binary 
companions, and then use the observed photometry to 
determine the posterior probability distribution of com¬ 
panion mass for each star in our sample. 

For simplicity, for each observed bulge star we consider 
11 discrete options J, J = 0,..., 10. The option J = 0 
implies a single star, J > 0 implies a binary system with 
mass ratio qj = J/ 10. The prior probability Prj of each 
option is consistent with the DM distribution with an 
overall binary fraction of 50%; thus Pro = 0.5, Prr_io = 
0.07, 0.072, 0.07, 0.063, 0.06, 0.05, 0.04, 0.035, 0.025, 
0.015. 

For each value of J, the total system mass ( Mt)j = 
( Mi)j + (A/ 2 ) j is chosen so as to match the total flux in 
the F&14W filter, using the appropriate mass-luminosity 
relation for age, metallicity and distance for both com¬ 
ponents (or only one component if J = 0). We then 
compute the likelihood Pd\m = P(Data|Model) of the 
measured total flux in the F814W filter, given the model, 
using a Gaussian distribution for the flux with the realis¬ 
tic photometric errors derived above. To the photometric 
errors derived from the AS tests, we added errors due to 
the presence of a metallicity spread, differential redden¬ 
ing and depth effects. These have been derived by using 
mass-luminosity relations for different metallicities, and 
by varying the distance modulus of 0.2 mag, and redden¬ 
ing of 0.1 mag, as described in §5. To each observed star 
we thus assign a probability distribution function (PDF) 



Fig. 8.— IMFs of the Galactic bulge derived by assuming differ¬ 
ent fraction of binaries and a DM mass-ratio distribution for the 
binaries. 


of the component masses over the allowed values of the 
mass ratio between components, according to the clas¬ 
sic Bayes formula: Pj = Prj * P D \ M / P(Data), where 
P(Data) is a normalization factor chosen to take into 
account the estimated completeness correction. 

Finally, we generate multiple realizations of the full 
stellar mass function by randomly drawing stellar dis¬ 
tributions with the probabilities thus determined. This 
procedure allows us to better understand and quantify 
the uncertainties arising from the correlated nature of 
the probabilities for each object (e.g., only one value of 
the mass ratio can be selected for each system). 

In practice, larger mass ratios qj generally correspond 
to redder P606W — P814VF colors at given F814TT-band 
flux; thus stars that lie red-ward of the main sequence 
of Fig. 4 will generally favor larger mass ratios, while 
stars located near the main sequence will be consistent 
with a single star or a low-mass binary companion which 
contributes little to the total flux. However, note that 
for many stars the photometric error is large enough that 
photometry (through the term Pd\m) does not provide a 
strong discriminant; for such cases, the final probability 
Pj for each option is the same as the prior probability 
Prj. By not taking into account the photometric color 
information, for instance, the distribution changes by ~ 
2-7 % in the VLM range, and by less than 1% at higher 
masses, i.e. M > 0.5 M 0 . 

As discussed in the following subsection, undetected 
binaries have a substantial impact in the inferred mass 
function, especially below « 0.5 M 0 . However, we must 
remark that the distribution of binary properties is un¬ 
certain and poorly constrained by the data at hand; 
changing the assumed binary fraction and the a priori 
distribution of mass ratios would also alter the derived 
mass function, as we show in Section 5.2. 
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The treatment above is somewhat simplified in com¬ 
parison with a fully Bayesian approach, in which we 
would consider fully the uncertainties in in the parame¬ 
ters of the model (metallicity, distance, reddening vari¬ 
ation), using for each an appropriate distribution rather 
than the “best” values. We defer this more complex and 
computationally expensive approach to the analysis of 
the full data set, including one more season of photome¬ 
try and eleven additional fields. 


5.2. Discussion 

One of the realizations of the Galactic bulge IMF is 
shown in Fig. 7. Error bars also include the uncertainties 
that come from statistical noise in the star counts. We 
generated 10,000 realizations of the same mass function 
and fitted them by adopting two power laws. The fit 
was performed by varying the mass break-point in the 
range —0.2 < log(M/M q ) < —0.3 and the lowest chi- 
square fit resulted with for a value of log{M/ Mq) = - 
0.25 (M = 0.56 Mq). The best estimate of the power- 
law slopes are a = —2.41±0.50 (dotted line) for higher 
masses, and a = — 1.25±0.19 (dashed) for lower masses, 
where dN/dM = CM a . 

We also fitted the Galactic bulge IMF by using a log¬ 
normal function described by the functional form: 


£(logM) = C exp 


[log{M) - log{M c )] 2 | 

2a 2 J 


( 2 ) 


with M c = 0.25±0.07 and a = 0.50±0.01 (solid green 
line in Fig. 7). 

The power-law slope for the high-mass range (M > 
0.56 Mq) agrees very well with the Salpeter IMF (a = 
—2.35) derived for solar neighborhood stars in the mass 
range 0.3 - 10 Mq (blue dashed-dotted line in Fig. 7). 
The log-normal mass function derived for disk st ars closer 
than 8 pc in the mass range 0.08 to 1.0 Mq bv lOhabrieH 
(I2003L(2005h has M c = 0.20 and a = 0.55 (dashed-double 
dotted red line) and agrees well with our Galactic bulge 
IMF. 

We derived the IMFs using the same method described 
in the previous section for different values of binary frac¬ 
tion, assuming a flat distribution of mass ratios, and a 
distribution given by DM. Fig. 8 shows one realiza¬ 
tion for each of the IMFs derived for different binary 
fractions and the DM mass-ratio distribution. In gen¬ 
eral, the IMF has two distinct slopes in the low- and 
high-mass ranges, and the slopes have only a weak de¬ 
pendence on the assumed mass-ratio distribution for the 
binaries. If we assume the DM mass-ratio distribution 
for the binary components, the slopes of the IMF at 
higher masses are —2.25, —2.36, —2.41 and —2.53, for a 
bulge binary fraction of 0, 30, 50 and 100%, respectively. 
If we assume a flat mass-ratio distribution, the slopes 
change only by 1 to 4% for binary mass fractions of 0 to 
100%. The corresponding slopes in the low-mass range 
are —0.89, —1.12, —1.25 and —1.51 for the DM mass-ratio 
distribution, and they change by 3 — 4% for a flat mass- 
ratio distribution. Full details including the error bars 
are given in Table 1. These results indicate that the 
effect of the presence of unresolved binaries is more pro¬ 
nounced in the low-mass range (~50%), than in the high- 
mass range (~12%). In the rest of the discussion, we use 
the IMF derived by assuming a binary fraction of 50% 


and a DM mass-ratio distribution. As discussed in §5.1, 
the presence of a substantial fraction of He-core WDs in 
the Galactic bulge suggests that the fraction of binaries 
in the bulge is larger than 30%. 

6. COMPARISON WITH OTHER INITIAL MASS 
FUNCTIONS 

6.1. Galactic bulge 

The Galactic bulge mass function was first measured 
by H098 based on a set of observations of the Baade’s 
window ( l = 1 °,b= —4°) collected with the Wide Field 
Planetary Camera 2 (WFPC2) on board HST. These 
data allowed them to derive a luminosity function down 
to F81AW ~ 24.3, corresponding to M ~ 0.3 Mq. No in¬ 
formation on proper motions was available, so disk stars 
are included in their study. But they applied a correction 
for the presence of unresolved binaries and found that the 
IMF of the bulge has a power-law slope of a = —2.2 in 
the high-mass range. The slope of the IMF flattens at 
~ 0.7 Mq, with a = —0.9 for a fraction of binaries of 
0% and —1.3 for 50% (see Table 2). H098 result for 
an assumed fraction of binaries of 50% agrees quite well, 
within the uncertainties, with what we obtained in our 
analysis for the same assumption on binaries, but the 
changing of power-law slope occurs at lower masses, ~ 
0.56 Mq, in our bulge IMF. 

A second study on the Galactic bulge mass function 
was published by ZO00, based on a set of observations 
collected in the F110W and F160W filters with NIC- 
MOS on board HST, covering a 22.5” x 22.5” field of 
view in a region of the bulge South of the Baade’s win¬ 
dow (l = 0.277°, b = —6.167°). To convert magnitudes to 
masses they used a mass-luminosity relation based on the 
same stellar models adopted in this investigation. They 
also did not have propel-motion information to separate 
bulge from disk stars, nor did they apply a correction 
for the presence of unresolved binaries. However, ZO00 
applied an overall reduction of the luminosity function 
by 11% for magnitudes brighter than J < 17, to take 
into account the contamination by disk stars. By fitting 
the IMF with a single power law they obtained a slope of 
a = —1.33±0.07, over the mass range 0.15 < M/Mq < 
1.0, while by using two different power laws they ob¬ 
tained a = —2.00±0.23 for masses M > 0.5 Mq, and 
a = —1.43±0.13 for lower masses (see Table 2). If we fit 
our IMF by using a single power law for the entire mass 
range (0.15 < M/Mq < 1.0), we obtain a range of slopes 
from a = —1.14±0.10 for no binaries to a = —1.56±0.10 
for 100% of binaries. The slope of the IMF not corrected 
for the presence of unresolved binaries is then shallower 
compared to the slope of the IMF obtained by ZO00 
(—1.14 vs —1.33). Moreover, the same IMF shows a 
much shallower slope in the low-mass regime compared 
to ZO00 mass function (—0.89 vs —1.43). This discrep¬ 
ancy might be due to the residual contamination by disk 
stars of ZO00 sample. Part of the difference could also be 
due to an intrinsic difference of stars observed by ZO00 
and stars in the SWEEPS field. From spectra collected 
by our group the stars in this region of the bulge have a 
similar metallicity distribution as the stars in the Baade’s 
Window, with main peaks at \M/H] ps —0.4,0 and 0. 3 
(|Hill et al.1 1201 H iBensbv et all 120131 : iNess et ah! I2013D . 
The metallicity distribution of the region of the Galactic 
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bulge observed by ZOOO shows a decrease in the fraction 
of metal-rich stars, with the average metallicity decreas¬ 
ing from [Fe/H\ ~ +0.03 in the Baade ’s window down 
to [Fe/H] ~ —0.12 (lZoccali et al.li200 &). However, such 
a small difference in the metallicity distribution cannot 
account for a ~ 20% difference in the IMF slope. 

6.2. Galactic disk 

The Galactic disk mass function has been constrained 
in the low-mass regime down to the hydrogen-burning 
limit and in_the brown dwarf regime by various studies. 
iSalpeteii ( 1955 ) derived the “original mass function” for 
solar neighborhood stars in the range 0.3 < M/Mq < 10 
and fitted it by using a single power law with a slope of 
a = —2.35. Later studies found that the disk mass func¬ 
tion can be reproduced either by a segmented power law 
or by a log-normal function. Table 2 lists the power law 
slopes and the characteristic masses and sigmas used by 
different stud i es to fit the G alactic d i sk ma ss function. 
iKroupa et al.l (I1993T) and later iKroupal (1200111 derived the 
IMF for disk stars within 5.2 pc of the Sun by taking 
into account a correction for the presence of unresolved 
binaries and fitting it with a power law with a slope of 
a = —2.2+0.3 in the mass range 0.5 < M/Mq < 1.0 
and of a = — 1.3+0.5 in the range 0.08 < M/Mq < 0.5. 
iGould et al.1 (119971) based their study of the Galactic disk 
mass function on photometry collected with the WFPC2 
and WFPC1 on board HST. They observed a sample of 
337 stars distributed in different regions of the disk and 
found a mass function with a slope close to Salpeter, 
a ~ —2.2, for masses in the ran ge 0 .6 < M / Mq < 1.0, 
and a ~ —0.9 for lower masses. IReid et all (l2002lf ob¬ 
served a sample of 558 main-sequence stars in the solar 
neighborhood in the mass range 0.1 < M/Mq < 3.0 and 
found that a power law with a slope of a = —1.3 fits the 
mass function in the lo w-ma ss range, i.e. for stars with M 
<0.7 Mq. iChabrieJ (1200511 adopted a log-normal func¬ 
tion to fit the Galactic disk IMF for single stars in the 
mass range 0.08 < M/Mq < 1.0, and found a character¬ 
istic mass M c = 0.20+0.02, and a =0.55+0.05. 

More recent analyses based on the Sloan Digital Sky 
Survey (SDSS) and the Two micron all sky survey 
(2MASS) data confirmed previous results, showing that 
the Galactic disk mass function can be reproduced ei¬ 
ther by assuming a segmented power law with slopes of 
a = —2.04/ —2.66 and a = —0.8/—0.98, for the high- and 
low-mass range, respectively, or by a l og-n ormal function 
with M r = 0.20/0.50 a = 0.22/0.37 (IGovev et al l 120081: 
IBochanski et al. 1201011 . 

The IMF we derived for the Galactic bulge is in 
very good agreement, wi thin uncertain ties, with the 
mass function obtained by IKroupal (|20011 f and iChabrieJ 
(2005) for the disk. On the ot her hand, the mass func¬ 
tions derived fo r the disk by IGovev et al.1 (12008H and 
IBochanski et al.l (I2010H have a slightly shallower slope 
compared to our bulge IMF in the low-mass regime (see 
Table 2), although the two mass functions would agree 
at lower masses by assuming the presence of no binaries 
in the bulge. 

7. THE STELLAR MASS-TO-LIGHT RATIO OF 
THE GALACTIC BULGE 

The stellar mass-to-light ratio ( M/L ) is an impor¬ 
tant parameter of a stellar population and depends on its 


TABLE 1 

Power-law slopes of the IMFs derived by assuming 

DIFFERENT BINARY FRACTIONS AND MASS-RATIO DISTRIBUTIONS 

for the Galactic bulge. 


Binary fraction 

Mass-ratio 

^Hiqh 

CX-Low 

0 


-2.25i0.50 

—0.89i0.20 

30 

DM 

—2.36i0.51 

— 1.12i0.19 

50 

DM 

—2.41i0.50 

— 1.25i0.19 

too 

DM 

—2.53i0.51 

—1.51i0.20 

30 

Flat 

—2.39i0.51 

—1.16i0.19 

50 

Flat 

—2.45i0.51 

—1.29i0.19 

too 

Flat 

-2.62i0.52 

— 1.55i0.20 


IMF. We used the mass function derived in this work and 
the total luminosity of stars observed in the SWEEPS 
field to estimate the stellar M/L of the Galactic bulge 
in the F814W and the F606IU filters. We obtain a to¬ 
tal mass for bulge stars in the SWEEPS field included 
in the mass range adopted to estimate the IMF, 0.16 
< M/Mq < 1.0, of 137,527+23,400 Mq. By extrapolat¬ 
ing the IMF with a power-law slope of a = —1.25 down 
to the hydrogen burning limit, we get an extra mass of 
14,310+2,400 Mq, for a total mass of « 152,000+23,500 
Mq included in the 0.10 < M/Mq < 1.0 mass range. 
Uncertainties take into account the error budget of the 
derived IMF. A constant mass of 1.0 Mq is assumed for 
bulge sub- and red-giant stars and red clump stars, for 
a total mass of 4,116 Mq. We do not take into account 
the mass loss along the RGB, but since the total mass 
of the giants is already very small compared to the mass 
of the MS stars, this has no effect on the final derivation 
of the mass-to-light ratio. We then assume that the IMF 
of the Galactic bulge has a constant Salpeter power-law 
slope for masses larger than 1.0 Mq and up to 120 Mq, 
and we integrate the IMF to obtain the number of stars 
that formed in this mass range. To estimate the mass 
currently in stell ar rem nants iii the b ulge we follow the 
prescriptions of iPercival et al.l (1200911 : stars with mass 
(i) 1 < M/Mq < 10, the remnant is a white dwarf; (ii) 
10 < M/Mq < 25, the remnant is a neutron star, and 
(iii) M > 25 Mq, the remnant is a black hole. In order 
to estimate the mass of white dwa rf re mna nts, we use d 
the initial-to-final mass relation bv lSalaris et al.l (120091 ). 
Mf = 0.084 Mi + 0.466 for initial masses less than 7 Mq 
and a constant final mass of 1.3 Mq for initial masses in 
the range 7 < M/Mq < 10, obtaining a total mass in 
white dwarfs of 53,912+9,200 Mq. For neutron stars we 
assumed a constant mass of 1.4 Mq and for black holes 
a mass equal to 1/3 of the initial mass, obtaining to¬ 
tal remnant masses of 3,905+600 and 11,151+1,900 Mq 
for neutron stars and black holes, respectively. By using 
the aforementioned values we found that the total stellar 
mass in the bulge SWEEPS field is M = 228,814 ±25,300 
Mq. 

We estimated the flux emitted by bulge stars in the 
SWEEPS field by using the proper-motion cleaned pho¬ 
tometric catalog corrected for the total fraction of stars 
and for completeness. We thus obtained a total lumi¬ 
nosity of Lp&uw ~ 104,000±2,000 Lq and Lpmew ~ 
71,000+1,400 Lq. 

The stellar mass-to-light values based on our IMF and 
the photometric catalog for the SWEEPS field are then 
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M/Lp saw = 2.2±0.3 and M/ Lfgo6W = 3.2±0.5. 

We estimated the stellar mass included in our field 
by also using two other assumptions for the mass dis¬ 
tribution at masses larger than 1 M 0 : constant slopes 
of a =-2.0 and of a =-2.7. In the first case, we obtain 
a larger total stellar mass, M = 254,505 ±28,100, and 
larger mass-to-light values, M/Lfsiaw = 2.4±0.4 and 
MlLF 606 W = 3.6±0.6, while in the second case we ob¬ 
tain smaller values, M = 219,079 ±24,200, M/Lpsuw = 
2.1±0.3 and M/Lp 606 W = 3.1±0.5. The total mass of 
the observed field and the stellar mass-to-light values 
estimated for the different cases are listed in Table 3. 

Finally, we also explored a more theoretical route 
and we computed the average bulge luminosity in the 
SWEEPS field by using two diffe rent synthetic popula¬ 
tion codes bv lCignoni et all (12013H and BASTI. For both 
simulations we generated a fake stellar population with 
properties resembling those in the Galactic bulge: so¬ 
lar metallicity, constant star formation rate between 12 
and 10 Gyr, our IMF, a binary fraction of 50%, distance 
modulus of 14.45 and reddening E(B — V) = 0.5. 

In the fir st case we used th e latest PARSEC stel¬ 
lar models (Bress an et al.ll20 12E Simulations were run 
until the number of stars in the magnitude range 20 
< F606W < 22 matched the observed number (~ 25000 
stars). This experiment was repeated 1,000 times. We 
found average values of Lpeo 6 W ~ 58,900 and Lpsi 4 W ~ 
92, 800. In order to evaluate the effect of metallicity dis¬ 
persion we also tested different Z values, namely 0.008 
and 0.03, corresponding to the metal-poor and metal- 
rich peaks of the metallicity distribution of the consid¬ 
ered field. In these cases we found Lp 606 W ~ 70,800 
and Lpsuw ~ 101,600 for the former metallicity, and 
Lfgosw ~ 53,800 and Lfsi 4 W ~ 87,500 for the lat¬ 
ter. As expected, lowering the metallicity causes an 
increase in the luminosity of the system. Interestingly 
enough, luminosity values estimated for the lower metal¬ 
licity, Z = 0.008, agree quite well with the observed val¬ 
ues, while values for the higher metallicities are system¬ 
atically lower than our flux estimates. A part of this dis¬ 
crepancy may be due to the possibility that a few very 
bright thin-disk stars are still contaminating our data, 
raising the inferred luminosities. In addition, the actual 
PARSEC models miss the asymptotic giant branch stel¬ 
lar phase, hence the predicted luminosities are likely to 
be underestimated. 

We repeated the same experiment using the BASTI 
models for the three different metallicites, and obtained 
Lf 606 W and Lpsi 4 W values as ~ 74,500 and ~ 103,000 
for Z=0.008, ~ 62, 500 and - 89, 700 for Z=0.02, and 
~ 57, 200 and ~ 86,200 for Z=0.03. In this case the lu¬ 
minosity estimates for the lower metallicity are also in 
very good agreement with the observed values, while the 
luminosities obtained for the solar and higher metallicites 
are systematically lower. On the other hand, the lumi¬ 
nosity estimates for the three metallicites based on the 
two different sets of models agree very well. 

Summarizing, we found a stellar mass-to-light ra¬ 
tios included in the range 2.1< M/Lpsuw < 2.4 and 
3.1< MILF 606 W < 3.6 according to different assump¬ 
tion on the slope of the IMF for masses larger than 
1 Mq. These are likely to be slightly lower estimates of 
the real stellar mass-light budget of bulge since a few 
bright disk stars might still be contaminating our lumi¬ 


nosity estimate. These values agree quite well, within 
the uncertainties, with the estimates provided by ZO00 
in the Johnson V filter, M/Ly ~ 3.4, by using their IMF 
with a single slope of a =-1.33 below 1 Mq, and by as¬ 
suming a constant Salpeter IMF for stars more massive 
than 1 Mq. 

8. MICROLENSING OPTICAL DEPTH 

Several thousand microlensing events have been de¬ 
tected to date towards the G alactic bulge, mainly by 
the OGLE (llldal ski et al.ll2015f) and MO A collaboration s 
(jBond et all 120011 : iSako et all 120081 : iSumi et all 120131 ). 
These microlensing events have been used by several in¬ 
vestigators to derive the total mass budget as well as the 
mass function of t h e lens es. 

iPaczvnski et all ( 1994 ). based on a small number of 
microlensing events, noticed that the observed microlens¬ 
ing optical depth is in excess of the theoretical esti¬ 
mates, indicating a much higher efficiency for microlens¬ 
ing by either bulge or disk lenses. This issue has been 
further investigated in re c ent years bv sev eral groups 
(IWvrzvkowski et ahl 12015 : Sumi et ahl 12013 ). A help¬ 
ful hint comes from the latest study by Wyrzykowski et 
al., which shows a dependence of the mean microlensing 
timescale on the Galactic latitude. This signals an in¬ 
creasing contribution from disk lenses closer to the plane 
relative to the height of the disk, which needs to be taken 
into account in the estimation of timescales and optical 
depths. 

Since the timescale of the microlensing event is pro¬ 
portional to the square root of the mass of the lens, 
the timescales can be used for a statistical estimate of 
the mass func t ion of the lenses. iZhao et all (119951) and 
lHan fc Gouldl (119961 ) used this approach and reported 
a mass function with a s lope of —2.0 and a cutoff at 
~ O.IMq. iCalchi Novati et ahl (120081) also attempted to 
fit the observed timescales of the microlensing events 
with a power-law distribution of lens masses and ob¬ 
tained a slope of —1.7 for the distribution. As pointed 
out by ZO00, there may be an extra bias in the observed 
timescales due to blending in the ground-based observa¬ 
tions, which causes the times scales to appear shorter 
than they actually are. This leads to an underestima¬ 
tion of the lens masses. Even so, the derived slope from 
microlensing observations is in between the two slopes of 
a = -2.41 and —1.25 derived here, and thus seems consis¬ 
tent. It would be interesting to extend this microlensing 
analysis to the currently available list of all the observed 
microlensing events. 

Finally, we note that the microlensing optical depth 
comes not only from the living main-sequence stars, but 
also from the white dwarfs, neutron stars and black holes. 
The mass-to-light ratio derived in this paper should help 
in deriving a more correct estimate of the microlensing 
optical depth. 

9. DISCUSSION AND CONCLUSIONS 

We have derived the IMF of the pure bulge compo¬ 
nent down to 0.15 M 0 . The Galactic bulge IMF can 
be fitted by two power laws, one with a steeper slope 
a = — 2.41±0.50 for M > 0.56 Mq, and another with a 
shallower slope a = —1.25±0.19 for the lower masses. A 
log-normal function fits the IMF too, with a characteris¬ 
tic mass of M c = 0.25±0.07 and cr = 0.50±0.01. 
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TABLE 2 

List of the different mass functions derived for the Galactic bulge and disk 


Reference 

Mass range 

a Hiah 

CX-Low 

Mbreak 

CX. 

M c 

a 

Galactic bulge 

This work 

0.15 - 1.0 

-2.41i0.50 

— 1.25i0.19 

0.56 


0.25i0.07 

0.50i0.01 

Holtzman et al. (1998) 

0.30 - 1.0 

-2.2 

-1.3 

0.7 




Zoccali et al. (2000) 

0.15 - 1.0 




— 1.33i0.07 



Galactic disk 

Salpeter (1955) 

0.30 - 10 




-2.35 



Kroupa et al. (1993, 2001) 

0.08 - 1.0 

-2.3F0.3 

— 1.3i0.5 

0.5 




Gould et al. (1997) 

0.08 - 1.0 

-2.2 

-0.9 

0.6 




Reid et al. (2002) 

0.10 - 3.0 




-1.3 



Chabrier (2005) 

0.10 - 1.0 





0.20i0.02 

0.55i0.05 

Covey et al. (2008) 

0.10 - 0.7 

-2.04 

-0.8 

0.32 

-1.1 

0.20 - 0.50 

0.22 - 0.37 

Bochanski et al. (2010) 

0.10 - 0.8 

-2.66i0.10 

—0.98i0.10 

0.32 


0.18i0.02 

0.34i0.05 


TABLE 3 

Stellar mass estimates and mass-to-light values for the 
Galactic bulge for different assumed IMF slopes for M > 

1 M 0 . 


a 

Stellar mass 

M/LpsiAW 

M/Lfqoqw 

Salpeter 

228,814 i25,300 

2.2i0.3 

3.2i0.5 

-2.0 

254,505 i28,100 

2.4i0.4 

3.6i0.6 

-2.7 

219,079 i24,200 

2.1i0.3 

3.1F0.5 


The slope of the IMF at high masses is mildly affected 
by the assumption on the fraction of unresolved binaries 
in the bidge or the distribution of their mass ratios. The 
high-mass slope ranges from a = —2.25 ± 0.50 for no 
binaries to a = —2.62 ± 0.52 for 100% of binaries in the 
bulge. On the other hand, the slope at lower masses 
changes significantly, ranging from a = —0.89 ± 0.20 for 
no binaries to a = —1.55 ± 0.20 for 100% of binaries. 

As we noted earlier, the slope of the IMF at the very 
low-mass range is crucial in estimating the mass budget 
of the Galactic bulge which contains ~20% of the mass of 
the Galaxy. Our deep HST observations obtained over 
a timescale of ~9 years allowed us to derive the mass 
function of the pure bulge component even in this low- 
mass range, which was previously not possible. 

The shape of the Galactic bidge IMF we derived in this 
work is in good agreement, within the uncertainties, with 
the IMFs derived previously by H098 for the Baade’s 
window, but our mass function extends to lower masses 
and it is purely based on bulge members with negligible 
contamination from disk stars. On the other hand, our 
IMF not corrected for the presence of unresolved binaries 
shows a slightly shallower slope compared to ZO00 IMF 
(—1.14 vs —1.33). This difference could be due to a small 
residual contamination by disk stars of the ZO00 sample, 
or to some intrinsic differences in the stars in the field 
observed by ZO00 and stars in the SWEEPS field. 

Our bulge IMF is in very good agreemen t wit h the mass 
func tion deriv e d for the Galactic disk bv iKroupal (120011 ) 
and iChabrierl (|2003f ) in the entire mass range, while it 
is steeper in the very low-mass regi me compared to the 
mas s functions derived for the disk bv lGould et, al.l (119971 ) 
and lReid et al.l ( 2002f). The P DMFs derived in the more 
recent studies of lCovev et al.1 (|2008l) and lBochanski et all 
(l2M0l) agree quite well with our IMF for the Galactic 
bulge in the high-mass range, but they still show a shal¬ 


lower slope in the low mass range. 

The characterization of the IMF in different stellar en¬ 
vironments is fundamental for investigating if the IMF 
has a dependence on th e ste llar metallicity and/or age. 
The recent work of iKalirai et alJ ( 201 3) showed that 
the IMF of the Small Magellanic Cloud (SMC, —1.5 < 
[Fe/H] > —1.0) is shallower than the Salpeter mass func¬ 
tion, a = —1.9, down to ~ 0.4 Mq, and does not show 
evidence for a tu r n-over in_the very low-mass regime. 
Furthermore, iGeha et al.l (120131 ) showed that the IMF 
of two metal-poor ([Fe/H] < —2.0) ultra faint galax¬ 
ies, Hercules and Leo IV, are even shallower, having a 
slope in the range a = —1.2 to 1.3 for masses larger 
than ss 0.5 Mq. In the higher-mass range (M > 0.5 Mq) 
where the mass function of these galaxies is well mea¬ 
sured, our bulge IMF is steeper than both the IMFs of the 
intermediate-metallicity environment of the SMC (—2.41 
vs —1.9) and the metal-poor environment of the ultra¬ 
faint galaxies (—2.41 vs —1.3 to —1.2), pointing towards 
a variation of the IMF with the global average metallicity 
of the stellar population. However, more data are needed 
to sample the IMF down to lower masses, i.e. 0.1 Mq, in 
the different environments, to confirm this preliminary 
result. 

We then used the derived IMF to estimate the stel¬ 
lar mass-to-light ratios of the Galactic bulge. We ob¬ 
tained for the two filters, values included in the range 
2.1 < M/Lfsiaw V 2.4 and 3.1 < M / Lfqoqw V 3.6 ac¬ 
cording to different assumption on the slope of the IMF 
for masses larger than 1 Mq. 

The shape of the mass function derived from microlens- 
ing observations has large uncertainties but is in consis¬ 
tent with the observed bulge IMF presented here. 
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